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We show that the ground state of the Bose fluid in a vessel-parallelepiped can be 
characterized not only by the quantum numbers p = (±2nl x /L x ,±2iily/Ly,±2nl z /L z ) 
flx, ly, h = 1, 2, 3, . . .) of an expansion in the eigenfunctions e tpr of a free particle, but also 
by the quantum numbers ki = (nl x / L x ,irl y / L y ,iil z / L z ) of an expansion in the eigenfunctions 
sin (h x x) sin (ki y y) sin (h z z) of a particle in the box. In the latter case, the one-particle 
condensate "is dispersed" over many lower levels with odd l x ,l y ,l z . 53% atoms of the 
condensate occupy the lowest level with k= (n/L x , 7t/L y , ir/L z ). The sum of condensates on 
all levels is determined by the same formula N c /N = Poo, as that for the condensate on the 
level with p = 0. The proposed k-representation supplements the traditional one (with the 
condensate on the level e = (hp) 2 /2m = 0) and allows one to consider the micro structure of 
a system from another side. 
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1 Introduction 

The problem of a condensate in the Bose fluid attracts attention as before. The condensate 
was predicted theoretically [T] — [3] and was then discovered experimentally. Moreover, the 
experimental values [H [5] agree sufficiently well with theoretical ones [3], E] |TT] . Intensively 
are studied the condensates in rarefied gases in traps. But some difficulties exist also. In 
particular, all the calculations of the condensate possess shortcomings decreasing the accuracy 
of results: sometimes, a lot of fitting parameters are used [Sj, or the correlation corrections are 
large [7J [10] ; in the most exact Monte Carlo simulations [6j [91 [11] , the interatomic potential 
is applied, though it is only effective at small distances and is known approximately. 

In addition, we meet the difficulty, which remains usually unnoticeable. A condensate in 
the Bose fluid is considered as the "macroscopic number of atoms with zero momentum." 
However, the fluid is placed always in a vessel, and the uncertainty relation hAp x Ax > h/2 
implies that the momentum of an atom cannot be zero due to the uncertainty Ap x > 1/2L. 
This point is obviously clear: the momentum is defined for translation-invariant systems, 
whereas the walls break the translational invariance, and the momentum of an atom is not 
conserved due to the collisions with walls. In addition, the momentum takes the discrete 
values p = 2n(j x / L x , j y / L y , j z /L z ) in the standard approach, and its uncertainty > 1/2L is 
comparable with the distance between adjacent values. Therefore, the collection of discrete 
levels is transformed in a collection of bands or even in a single band. This means that the 
momentum is a "bad" quantum number. However, this circumstance creates no significant 
problem, because we can consider p as an eigenvalue of the Hamiltonian (e = (hp) 2 /2m) and 



can use the traditional formalism. But it is possible to go along a nontraditional way with the 
construction of another approach to the description of the microstructure of systems. Such 
new approach will be considered in what follows. 



2 Structure of a condensate for the Bose fluid in a vessel 

Consider the Bose fluid in a vessel at T = 0. In the standard approach with periodic boundary 
conditions (BCs), the condensate n = N p=0 /N is determined from the formula 

_/V f e -ip(ri-r 2 ) 

/ dridr 2 p(ri,r 2 ) — 2 , (1) 

where N is the total number of identical Bose particles, V is the volume of the system, 
p(r l5 r 2 ) is the one-particle density matrix, 

p(r a , r b ) = V J dr 2 ... dr N %(r a , r 2 , . . . , r N )^ (r b , r 2 , . . . , r N ), (2) 

and \l/o is the wave function of the ground state of the system. For an infinite system, we 
have p(ri,r 2 ) = p(\ri — r 2 |) due to the translational invariance and the isotropy. Therefore, 

H ° = ~~/V^ = V* j dTldT *P(\ Tl ~ r2 D = P(\ Tl ~ r2 ' ^ °°) = Poo- (3) 
The presented definition of the condensate follows from the expansion 

#(r) = ( (ipa(p)e ipr (27r)- 3 (4) 



in the eigenfunctions e ipr of the momentum operator, which form the complete collection of 
orthonormalized basis functions of the continuous spectrum. The exponential function e ipr 
is an eigenfunction of the Hamiltonian and the momentum operator for a free particle in the 
infinite or finite closed volume. The closed ring systems are described by periodic BCs and 
are realized in the nature only in the one- or two-dimensional case. 

Consider the finite-size systems with nonperiodic BCs. In particular, the boundary con- 
ditions for the three-dimensional systems are always nonperiodic. As was noted in [12] , the 
expansion for the system in a box should be performed in an eigenfunctions of this system. 
Work [12J gives the solution for iV point Bose particles in a box-cube L x L x L character- 
ized by an interatomic potential in the form of the 5-function. The BCs were taken to be 
nV^ = 0. It was found that all atoms are on the lowest level at a repulsive potential. At 
an attractive potential, the atoms are distributed over (~ N 2 ^) lower levels, but with the 
microscopic (~ iV 3 / 5 ) population of each level. In this case, the condensate exists in the 
"generalized" sense [12] : the total number of atoms on those levels turns out macroscopic. 

But the real atoms are not point-like and possess a smoother potential, than the 5-function. 
Therefore, the wave functions may have a different structure as compared with that of the 
solutions in [T2]. We may expect also a different solution for the condensate. We now consider 
such realistic system of Bose particles in the box L x x L y x L z . For simplicity, let us consider 
firstly the ID case. We model the potential of a wall as a step with finite height: 



U w (x) 



U > x < 0, x > L , . 

< x < L. W 



Let the potentials of both walls be identical (their distinction will not lead to basic changes). 
The eigenfunctions of a free particle in such a box are known |14j : 



ci sin (5i)e K ' x , x < 

C\ sin [k\x + 5i), < x < L x (6) 

c l (-l) l+1 sin (5i)e- Kl ^ x - L \ x > L X} 



h = ^2mE l /h = vrZ/L, - 25t/L x , (7) 
ki = ^2m(U - Ei)/h > 0, 8 t = arcsin (7^/^/2), (8) 



Here, Z = 1,2,3, .. ., 7 X = 77- y an d the values of arcsin are taken between and n/2. 

The boundary conditions involve the continuity of ^i(x), ^[(x) and are more physical than 
the BCs nV^ = in [12]. For great systems, 7 X <C 1: for He 4 atoms at U = 100 K and 
L x = 10 cm, we have 7^ ~ 4.9 x 10~ 10 . For I <C l/7rj x , we obtain 

ki w Zfci = 1(1 - -f x )n/L x , 61 w Wi = /vr 7x /2. (10) 

Solutions © describe bound states with discrete spectrum. Here, we omit the solutions 
with continuous spectrum, which correspond to the infinite motion with E > U. The 
solutions for the continuous and discrete spectra are eigenfunctions of the Hermitian operator 
(Hamiltonian). By the theorems of quantum mechanics, these functions form the complete 
orthonormalized collection of basis functions, in which any function ty(x) with the appropriate 
smooth properties and given on the interval x e] — 00, oof can be expanded. Therefore, any 
solution for a free particle localized in a ID box (jSJ) can be written in the form 

00 

tf(aO = £}&jtf f (a:), h= I dx*{x)^{x), (11) 
' -00 

and the probability for the particle in the state ty(x) to have the wave number ki is equal to 
l^l 2 - 

Let us consider the Bose fluid in a vessel. In the presence of a vessel, the momentum 
is not a quite good quantum number. Therefore, we are not based on the momentum op- 
erator, while choosing the eigenvalues for the expansion of the density matrix. Then the 
exponents e tPxX lose the status of exceptional functions: now, we can use only the eigenfunc- 
tions of the Hamiltonian of a free particle, which are the collections of all possible linear 
combinations of the form ae tPxX + he~ tPxX with different p x , where p x are multiple to 2tt/L. 
Among those collections, we separate two ones: first, the exponential functions e iPxX (because 
they give solutions for the quasiparticles in the form of traveling waves, and the phonons 
in real systems are namely traveling wave packets, rather than standing waves; we call it 
p-representation), and, second, the sines sin (fyx + 5i) flUD (since they take the boundary con- 
ditions into account, by ensuring the decrease of the wave functions near the boundary; we 
call it the k- representation). Here, the wave vectors k and p are similar to the momentum, 
but they are not the momentum, but a quantum numbers of the Hamiltonian. None of these 
collections of eigenfunctions can be called definitely the best. In our view, the exponential 
functions describe better the quasiparticles, whereas the sines do the ground state and the 
condensate. 



In the p-representation, we obtain obviously the well-known results with the condensate 
on a level with p = 0. Of a significant interest is the k- representation based on the functions 
©• We now consider it in more details. The wave vector takes the values 

k = k l= ( nlx ~ 2Sl * nl y ~ 26l v nl * z 2Sl * \ (12) 



with l x ,l y ,l z — 1,2,3,... We omit the negative components, since the wave function of a 
particle for them differs only by a sign (these are equivalent states). For the Bose fluid states 
of the discrete spectrum, the number of particles with the wave vector k is 

oo 

^ = i J dridrapCrx.raJ^Cxx)*,.^)^^)*,.^)^^)*,.^). (13) 

— oo 

Let us verify the correctness of the coefficient V -1 . The property of eigenfunctions 

J2 *r(ri)*i(r 2 ) + f df^(ri)tf f (r 2 ) = 5{v x - r 2 ) (14) 
l 



and the normalization of \I/o yield 



k 



i.e., the coefficient is proper. The integral over the continuous spectrum is probably much 
less than 1. 

Let us determine N^. According to ({TBI) , we need to integrate over the whole volume, 
inside and outside of the box. But we do not know the solution for \l/ outside of the box. 
However, it is clear that this solution decreases exponentially, with the penetration distance 
~ 1 A for realistic walls. Therefore, we neglect the region outside of the box. Then, in the 



3D case, we obtain 



N y 



uXi'-'Vi'-'Z 

X - J dT 1 dr 2 p(Ti, v 2 )c 2 lx c 2 ly c 2 lz sin (ki x x x + 8 lx ) sin (k lx x 2 + 6 ix 



x sin (ki y yi + 5 iy ) sin {h y y 2 + 5 iy ) sin {h z zx + 5 t J sin (k iz z 2 + 5i x ). (16) 

We now have p(ri,r 2 ) 7^ p(ri — r 2 ), generally speaking. At distances from the wall much 
larger than the interatomic distance R, the properties of \&o are the same [13J as those for 
an infinite system. Therefore, at I - ! and r 2 far from the walls and |r x — r 2 | 3> R, the relation 
p(ri, r 2 ) = poo = const is true. With regard for it, relation (fl6|) yields 

_ ( Sc lx ci y c h cos 5 tx cos 5 ly cos 5 h \ 2 Poo 



n V Kh v K J v 



if l x ,l y ,l z are odd, and N kl /N = 0, if at least one of l x ,l y ,l z is even. Let us consider odd 
l x ,l y ,l z . The value of is large at small ki, when relations fflOl) hold. In this case, 

iVk, 8 3 poo 



Since ~ 1, we obtain the condensate ~ N for all levels with small l x , l y , l z . The sum 
of condensates on all levels 



EA k ^ 8 3 poo \ - 1 
k N~ ** . (1 + 2^(1 + 2^(1 + 2^ (i9j 



>0 



{jx,j y ,jz = 0, 1,2,3, . . .) is exactly equal to the condensate on the level p = for the p- 

representation. This interesting property is general. We denote V^ 1 ^ 2 = |r), p^ = p, 

^ i x (x)ty i y (y)ty i z (z) = |1) and expand the function \I/(r) , equal to V~ 1 / 2 in the box and to 

zero outside of it, in the functions |1) and the continuous-spectrum functions of a particle in 

the box. The latter enter the expansion with zero coefficients. Therefore, \I/(r) = c i 1 1} ^ 

l 

C\ = (r|l). Then, in the approximation p(rx,r 2 ) = p^, we have 

J2N p /N = N p=0 /N = Poo = I dr Poo (e ipr /V)\ p _ 
p J 

= (r|p|r) = (r|l)(l|p|m)(m|r) = £ ei<4(l|p|m) 

l,m 

= Ei c 'i 2 ( 1 i/ ) i 1 ) = ^Ei c ii 2 = E iV ^/ iV - (2°) 

1 1 k[ 

We do not use the specific form of the functions Therefore, the summary condensate 
N^JN is the same for any basis of eigenfunctions. This property is related to the fact that 

ki 

the constant p^ can be expanded in various bases. 

The real value of n c is somewhat less than that in (fl9l) . The condensate corresponds to 
the macroscopic filling of a level, but sum (fT9~!) includes also large values of j x ,j y ,jz, at which 
the filling is small. We call the level macroscopically filled, if > q c N. The value of q c is 
conditional, but the estimate q c ~ (1 — 100) x iV~ 1//3 seems reasonable. The numerical analysis 
indicates that, at q c = 0.001, we have n c « 0.914POO and the number of condensate levels 
iV c = 49. At q c = 10" 6 , we obtain n c « 0.994poo and N c = 3920. Whereas n c « 0.999poo and 
N c « 61400 at q c = 10~ 8 . For great systems iV 1/3 > 10 s . In this case q c ~ 10" 6 -10~ 8 , and n c 
differs slightly from n c = p^ (|T9l . 

It is clear that, for the ID and 2D cases, the properties of the condensate are analogous. 

For a different shape of the vessel, we need to use different basis functions at the expansion, 
but the summary condensate is not changed, according to fl2"U|) . If we consider the realistic 
potential of a wall (smoother than that of a step and like the interaction potential of two 
helium atoms), the summary condensate is not changed as well, but, probably, will be more 
strongly dispersed over levels. 

The additional information is given by the exact wave function \l/o [13] of a system of iV 
interacting Bose particles in a box-parallelepiped with impermeable walls (U = oo): 

= ^e 5 -' Y[ {sin (k lx xj) sin (k ly y,) sin (k lz z 3 )} , (21) 
j 

q+qi^o (i) 

S« = £^(q)p_ q+ £ ^i|^p qi p_ qi _ q + ..., (22) 

q^0 q, qi ^0 V7V 

N 

where ki = (k lx ,k ly ,k lz ) = {ir/L x , Tr/L y , n/L z ), p k = ^ e~ tVx \ and is the well-known 

i=i 

solution (15] for \l/ i n the thermodynamic limit: 

ln ^o =2^ ~^r~ PuP-n + 2^ 3! -7ft p^p^p-^-^ + ■■■ (23) 



The components of q and q^ are multiple to 2tt/L. If we expand both exponential functions 
in ( 121]) in a series, represent pk in terms of cosines and sines, and convolve them with sines 



from the bare product of sines in (1211) . we obtain the sum of products of N sines and cosines 
with wave vectors multiple to tt/L. This yields two conclusions for particles in a box. 1) 
The wave vectors of interacting particles "are quantized" in the same way like those of free 
particles. 2) If the interatomic interaction is switched-off, \l/o is reduced to a product of sines, 
i.e., to the solution for the ground state of N free particles, when all particles are in the 
condensate with k = ki. The structure of condensate (fT9l) agrees with this solution: if the 
interaction tends to zero, the atoms from above-condensate levels "fall" on the condensate 
levels; if the interaction is switched-off completely, the condensate levels become empty by 
jump: the atoms transit to the lowest level with k = ki. But if the condensate is defined in 
the standard way (jTJ), we cannot obtain iV particles on the level with k = kx at the switching- 
off of the interaction. In this relation, the standard definition of a condensate is worse. The 
condensate measured in experiment is, in fact, the quantity p^. 

Notice that the composite condensate (fT9|) belongs to the generalized condensates, which 
was considered in several works, but only for free particles or point-like potentials or cyclic 
BCs (see, f.e., the recent paper [T6]). 

Let us compare the formulas for the number of above-condensate atoms, N p , in the k- and 
p-representations for k,p 3> tt/L. Let us expand the density matrix in a Fourier series: 

p(r 1 -r 2 ) = i 7 ^/(q)e^ ri - r2 ). (24) 
q 

For the p-representation, we use periodic BCs. Then V* = V, and q runs the values 
2n(j x /L x ,j y /L y ,j z /L z ). Relation (P yields 

N p = f(p)N/V. (25) 

For the k-representation, V* = 8V, and q runs the values n(j x /L x , j y / L y , j z / L z ), and rela- 
tions ([HD and (J25) yield 

N k = nf(k)N/V, n = r lx r ly r lz , (26) 

1 cos 2 5i sin 2 Si 

ri = - H — H 

x 6 3 25? 

l x = 1, 2, . . . , /™ ax — 1. The value of l"^ ax = l/8\ can be determined from the condition 
Simax = 7r/2 and gives the maximum value of ki x . The analogous result holds for the y- 
and ^-components. The vectors k > ]& max correspond to the continuous spectrum. At 
k < k max /2 (for all components), r\ ~ 1 is true. At larger k, r\ decreases, as k increases, 
and ri (k max ) « 0.05. For He 4 atoms at U = 25 K, we have k max = VOmU/h w 7.2 A" 1 , 

which equals to the threshold momentum k m —fa 3.6 A of the dispersion curve. We do not 
exclude that this property explain the threshold. Thus, r\ ~ 1 for all k for which /(k) is not 
too small, and the formulas for N k and A^ p are identical. In view of the positiveness of the 
components of k, we have Yl = ^ for narrow bands Ak and Ap at p = k. Therefore, for 

k p 

k, p < k. max /2, the numbers of above-condensate atoms in the k- and p-representations are 
the same (£) N k w /(k)A^fc 2 AA;/(27r 2 )). 

k 

According to [13] . the solution for the wave function \l/o under the zero BCs differs from 
that under the cyclic BCs. In view of this, the values of must also be different. Hence, the 
amount of the condensate and the number of above-condensate atoms under the zero BCs 
must differ from those under the cyclic BCs, though the formulas are identical. However, we 
have not calculated in the present work. 



3 Conclusion 



The analysis shows that the presence of boundaries allow one to propose the new definition 
of the condensate. At this definition, the condensate is dispersed over many lowest atomic 
levels. This effect is related to that the walls of a vessel play the role of a resonator and affect 
the spectrum of eigenmodes of the system. 

The proposed k-representation is equivalent, on the whole, to the traditional ones, but 
it allows one to observe new properties of a system, which are hidden, if the ordinary p- 
representation is used. 
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